Dispersion of waves in relativistic plasmas with isotropic particle distributions 
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The dispersion laws of Langmuir and transverse waves are calculated in the relativistic non- 
magnetized formalism for several isotropic particle distributions: thermal, power-law, relativistic 
Lorentzian k, and hybrid j3. For Langmuir waves the parameters of superluminal undamped, sub- 
luminal damped principal and higher modes are determined for a range of distribution parameters. 
The undamped and principal damped modes are found to match smoothly. Principal damped and 
second damped modes are found not to match smoothly. The presence of maximum wavenumber is 
discovered above that no longitudinal modes formally exist. The higher damped modes are discov- 
ered to be qualitatively different for thermal and certain non-thermal distributions. Consistently 
with the known results, the Landau damping is calculated to be stronger for non-thermal power- 
law-like distributions. The dispersion law is obtained for the single undamped transverse mode. 
The analytic results for the simplest distributions are provided. 



I. INTRODUCTION 

The non-thermal distributions of electrons are as im- 
portant as thermal for astrophysical plasmas. Shocked 
and turbulent medium is likely to accelerate the elec- 
trons into the power-laws. Gamma- Ray burstsA (GRBs), 
jets from compact sources^, low luminosity active galac- 
tic nuclei 3 (LLAGNs) show evidence of the non-thermal 
relativistic distributions. The electrons are non-thermal 
and mildly relativistic near Earth4 and in solar corona 5 . 

The properties of EM waves propagating through such 
medium are worth knowing even in a non-magnetized 
case. EM waves can be generated in turbulence, propa- 
gate some distance and dump via the Landau damping. 
Thus energy is redistributed. Whereas the realistic tur- 
bulence is non-linear, only the linear waves are considered 
in this work. No back-reaction of waves on the electron 
distributions is assumed. 

The reasonable relativistic electron distributions are 
thermal (|10p. power-law ||TTJ), Lorentzian k (|12p and hy- 
brid (3 (I13[) . They are taken to be normalized to unity 
when integrated over momenta. I take all the quantities 
in the paper to be dimensionless for the sake of clarity 
and brevity. Temperature T of thermal distribution is 
evaluated in the units of mc 2 //cs for the particles with 
mass m. Here c is the speed of light and ks is Boltzmann 
constant. The inverse temperature 



P 



1 

T 



(1) 



is sometimes denoted as p in the literature. Momenta of 
particles p is measured in mc. The speed of light c is set 
to unity. Thus 
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= v / T 



p 



(2) 



is the particles dimensionless energy. The non-relativistic 
plasma frequency in CGS units (for particles with charge 
q) 



is employed to normalize the absolute values of a 
wavenumber k and frequency u> as£ 



K = — . 

UJ p 



(4) 



The scalar K is the absolute value of the vector K. 



The general characteristics of relativistic plasma are 
the following. Transverse waves have the phase veloc- 
ity greater than the speed of light fl/K > 1 and so are 
undamped. They only have the single mode. The longi- 
tudinal waves exhibit the wider range of phenomena^. 
The phase velocity goes from > 1 to < 1 as K increases, 
the undamped mode at small K matches the principal 
damped mode that exists at higher K up to K max & Be- 
sides the principal mode, the finite set of higher damped 
modes exists. The present paper elaborates on all these 
modes for various isotropic particle distribution. Either 
a mixture of electrons and positrons with same distri- 
butions or immobile ions are considered, so that the ion 
sound does not appear. 

The paper is organized as follows. The formalism of 
linear plasma dispersion is reviewed in fJTTJ The results 
of numeric evaluations are presented in mill Some ana- 
lytic formulas can be found in fJTV] I conclude with the 
discussion and the summary in SjVj 



II. LINEAR DISPERSION LAWS 



Waves propagating in plasma feel the plasma response, 
which can be characterized by the permittivity tensor 
e(K, Q) in the linear regime. Permittivity tensor^ 
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e(K,fi) = 1 



E 

species 



l 

«2 



15/ 



PP 



p dp 7 — K • p/f) 



d 3 p (5) 



2 



leads&iS after the integration over the polar angles to 
longitudinal permittivity 
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and transverse permittivity 
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where 70 = (1 — ti 2 / K 2 ) 1 / 2 . The Landau rule is applied, 
so that a = for Re[ft 2 ] > X 2 and for Re[f7 2 ] < K 2 one 
has£ 



for Im[n] > 0, 
2 for Im[fi] < 0. 



(8) 



The case with Im[fi] = for Rc[0 2 ] < K 2 is unphysical 11 . 
Only one sort of species is considered in the above for- 
mulas without the loss of generality. The plasma density 
n is that of the mobile species in the the case of immobile 
ions or the total density for electron-positron plasma. 



III. NUMERIC RESULTS 

The dispersion laws Q(K) of waves propagating in 
plasma are determined as solutions of the equations 
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Thermal distribution 



The thermal relativistic distribution of particles 



Mp) = 



pexp(-p7) 
4ttK 2 (p) 



(10) 



describes highly collisional relaxed plasma. Here and be- 
low K n represents n-th modified Bessel function of the 
second kind. 

The longitudinal waves in thermal plasma have variety 
of features. The non-relativistic theory-— predicts the in- 
finite number of damped modes, whereas the relativistic 
analysis 7 limits the damped modes to a finite set. In ad- 
dition, the superluminal undamped part of the spectrum 
appears. The boundary K for these modes are shown on 
Fig. [1] These are the solutions of the dispersion relation 
@ for real f2 and K. The transition between the damped 
(high K) and undamped (low K) modes is indicated by 
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FIG. 1: (Color online) Upper solid (red) - maximum K for 
principal damped longitudinal harmonic, lower solid (blue) 
- for undamped harmonic; dashed (green) - K boundaries 
for second damped harmonic, dot-dashed (cyan) - for third 
damped harmonic. 



the lower solid blue line. The principal damped mode 
exists for K limited by the upper solid red line, whose K 
grows exponentially fast as T — + 0. No modes formally 
exist above that line, as the ion sound is absent. How- 
ever, the calculations in this large- K area are unphysical, 
because the set of particles constituting plasma does not 
behave coherently at very low wavelengths, as the wave- 
length approaches the particle effective mean free path 
about Debye radius Xd- Thus the branch ceases to exist 
at wavenumbers k > 1/Au 13 . Every realization of the 
disribution function f(p) gives rise to different evolution 
of the imposed initial condition. The dashed green lines 
indicate the upper and lower liming K for the second 
damped mode. This mode exists for any temperature. 
Its region of allowed K overlaps with that for the princi- 
pal mode. The lower limiting K goes to at p ~ 2 and is 
zero for higher p (lower temperature). The same is true 
for the third damping mode (boundaries in dot-dashed 
cyan): the lower limiting K goes to at p ~ 8. The third 
damped mode exists only for p > 6.7. 

The mode completion effect was claimed to exist by 
Schlikeiser— . It consists of the smooth transition of 
Re(f2(.ftr)) between the principal and second damped 
modes. This is quite surprising given the fact that 
Im(Q(K)) is discontinuous, thus Q(K) is discontinuous. 
The careful analysis shows that the real part is also dis- 
continuous. The real and imaginary parts of n(K) for 
p = 2 are shown on Fig. [5] (see Fig. 3 in Ref. (0) for com- 
parison). In agreement with ^-overlapping of the princi- 
pal and second damped modes (Fig.[T]), K max w 1.060 for 
the second damped modes is larger than K m i n w 1.000 for 
the principal damped mode. Thus, the mode completion 
effect does not occur. 

The isocontours of constant Re(fi) and — Im(f2) for 
the principal mode are shown on, respectively, Fig. [3] 
and Fig. 2J These are the extensions of Figs. 6 and 
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FIG. 2: (Color online) Dispersion relation for longitudinal 
waves for p = 2. Undamped mode is dot-dashed (blue), princi- 
pal damped mode is dark solid (black) , second damped mode 
is light solid (cyan). 

7 in Ref. (@) to higher temperatures (lower p) with 
the maximum K boundary (see Fig. [T]) employed. No 
modes formally exist below the dark solid line denoted 
by "max(K)". The modes are undamped to the left from 
the light solid cyan line "fl — K, Im(f2) = 0" and damped 
to the right. The plasma frequency at the infinite wave- 
length is shown in dashed green. The correspondent plot 
for the transverse waves is given in Ref. (6) along with 
the approximations. 

B. Power- law distribution 
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FIG. 3: (Color online) Dispersion contours for longitudinal 
waves, showing Re(f2) in (K, log 10 p) plane (solid blue lines), 
zero damping boundary (light solid cyan line), maximum K 
curve (dark solid black line), Q as a function of p for K — 
(dashed green line). 



Isocontours of constant -Im(fi) 
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FIG. 4: (Color online) Dispersion contours for longitudinal 
waves, showing — Im(f2) in (K, log 10 p) plane (solid blue lines), 
zero damping boundary (light solid cyan line), maximum K 
curve (dark solid black line), Q as a function of p for K — 
(dashed green line). 



The distribution 



fp(p) 



T(k) 



7r 3 / 2 r( K - 1) 



(ii) 



represents a simplest f(p). Also, it produces many ana- 
lytic results (see PV[). It is often used to interpret the 
results of astronomical observations of jets^ or any object 
where shock acceleration takes place. In applications this 
distribution is usually applied with the limited range of 
7, but for my calculation the entire range of 7 is taken. 

The migration of critical points of longitudinal disper- 
sion relation is shown on Fig. The superluminal mode 



smoothly converts to a damped principal mode at crit- 
ical K about 1, shown in solid blue on the figure. The 
maximum K of the principal mode (dashed black curve) 
grows exponentially fast as a function of k similarly to the 
behavior of maximum K boundary for the thermal distri- 
bution Fig. [T] The behavior of higher damped modes is, 
in contrast, different. The second damped mode appears 
at k — 4.58. At higher n second damped mode can as- 
sume K from to a critical K, shown in dot-dashed green 
on the figure. The presence of the second and higher or- 
der damped modes at K = is a feature of the power-law 
distribution and is not observed for the thermal distribu- 
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FIG. 5: (Color online) Maximum K for principal damped lon- 
gitudinal mode (dashed black line), for second damped mode 
(dot-dashed green line), for undamped mode (solid blue line). 

The isocontours of constant Re(f2) and — Im(f2) for 
the principal mode are shown on, respectively, Fig. [6] 
and Fig. [71 The choosen small range of k <G [2,3] re- 
flects the astronomical need to consider distributions 
N(j) ^ j- x d-y for 7 > 1 with x € [2,4]. As in the 
case of thermal plasma, no modes formally exist below 
the solid black line denoted by "max(K)". The modes 
are undamped to the left from the light solid cyan line 
v fl = K,Im(r2) = 0" and damped to the right. The 
plasma frequency at the infinite wavelength is shown 
in dashed green. The isocontour K appears to depend 
only weakly on k, however the isocontour — Im(f2) val- 
ues (Fig.[jj are at least 1.5 times larger for the same K, 
than in case of thermal distribution (Fig. [4|) . The com- 
parison of damping coefficients for thermal and hybrid 
distributions is given in the next subsection. 

C. Hybrid distribution 

I choose the relativistic isotropic Lorentzian k 
distributional 

/«(p)oc(l + ^^)^ +1) , (12) 

which tends to relativistic thermal distribution with 
T = 2 as k approaches infinity. The normalization 
coefficient is calculated numerically from the condition 
J °° 4:Trp 2 fn(p)dp — 1. This distribution is the realistic 
choice when the acceleration of particles is balanced by 
the relaxation processes: almost thermal distribution at 
high k and almost power-law at low k. Thus, it is possi- 
ble to compare the wave propagation for thermal, slightly 
non-thermal and highly non-thermal distributions. 

The comparison of distributions is made on Fig [5] The 
inverse temperature p for fx(p) and the temperature- 
like parameter 2 for fu{p) are choosen in such a way 



Isocontours of constant Re(fi) 




FIG. 6: (Color online) Dispersion contours for longitudinal 
waves, showing Re(f2) in the {K,k) plane (solid blue lines), 
zero damping boundary (light solid cyan line), maximum K 
curve (dark solid black line), O as a function k for K = 
(dashed green line). 



Isocontours of constant -Im(fl) 




FIG. 7: (Color online) Dispersion contours for longitudinal 
waves, showing — Im(f2) in the (K,k) plane (solid blue lines), 
zero damping boundary (light solid cyan line), maximum K 
curve (dark solid black line), f2 as a function k for K = 
(dashed green line). 



that the average kinetic energy of particle Ek — J (7 — 
l)4irp 2 f(p)dp is equal to 1. The choice of constant Ek = 
1 is made to show the typical behavior of dispersion re- 
lations for plasmas with approximately the same ener- 
getics. The thermal distribution is shown in dot-dashed 
black. The transition from light green to dark blue 
corresponds to the increase of k through the discrete 
set k = 4.5, 5, 5.5, 6, 7, 10, 20. The relativistic Lorentzian 
k distribution appears to have more particles at lower 
7 w 1.2 and higher 7^5, but fewer particles in the in- 
termediate region 7 ~ 3. For the contrast I added the 
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Relativistic Lorentzian *:, thermal and hybrid ji distributions 




FIG. 8: (Color online) Distributions of particles with kinetic 
energy Ek = 1 : thermal [n — > oo) (dot-dashed black line), 
Lorentzian ft with k = 4.5, 5, 5.5, 6, 7, 10, 20 (light green to 
dark blue lines), hybrid (3 distribution (infinite Ek) (dashed 
red line). 

hybrid distribution 

W) = ?(iTW (13) 

with (3 — 2 that is shown by the dashed red line. It was 
not normalized to Ek = 1, since it has the divergent 
total energy. However, the dispersion laws exist for it. 



Re(fl) for longitudinal waves 
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FIG. 9: (Color online) Re(£l(K)) dependence for longitudinal 
waves for Lorentzian ft distribution: limiting thermal case 
(k — > oo) (dot-dashed black line), various ft from 4.5 to 20 
(light green to dark blue lines), hybrid (3 distribution (dashed 
red line). 

The dependencies of Ke(il(K)) for longitudinal waves, 
— lm(il(K)) for longitudinal waves and fl(K) for trans- 
verse waves are shown for principal modes on, respec- 
tively, Fig. [H Fig. [TUl Fig. [TTJ The thin black line on 
Fig. [5] separates the undamped mode (on the left) from 



-Im(fl) for longitudinal waves 



1.0 




1.0 1.5 2.0 2.5 3.0 

K 



FIG. 10: (Color online) — Im(Q(K)) dependence for longitudi- 
nal waves for Lorentzian ft distribution: limiting thermal case 
(ft — > oo) (dot-dashed black line), various ft from 4.5 to 20 
(light green to dark blue lines), hybrid (3 distribution (dashed 
red line). 



fl for transverse waves 




K 

FIG. 11: (Color online) Ke(Q(K)) dependence for transverse 
waves for Lorentzian ft distribution: limiting thermal case 
(k — + oo) (dot-dashed black line), various ft from 4.5 to 20 
(light green to dark blue lines), hybrid (3 distribution (dashed 
red line). 

the damped one (on the right). The long- wavelength be- 
havior of longitudinal modes is predictable: more mobile 
lower 7 species in thermal distribution have the lowest 
plasma frequency. By the same reason f2 at K = for 
transverse waves (Fig. II 1 [) is larger for smaller n. The 
variation of Re(S!) is within 20% for both longitudinal 
and transverse waves. The Landau damping (Fig. HOj) 
shows larger variation of up to 3 times. The electrons 
responsible for Landau damping mainly have high 7^5, 
what makes distributions with smaller k dissipate waves 
more effectively. The similar result is observed for non- 
relativistic Lorentzian n distribution^. The behavior of 
plasma waves for hybrid [3 distribution is very different 
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from that for thermal or Lorentzian k distribution. 



General case 



IV. ANALYTIC FORMULAS 



The permittivity components are 



Thermal distribution 



The Trubnikov's form 1 - of the dielectric tensor reads 



where 



da 



K 2 {r) 



K 2 



Suv — Oi 



2 K^Ky K 3 (r) 



n 2 



p 2 — Hap + a 2 ( — 1 



-,1/2 



I denote 



1 



n 2 



(14) 
(15) 

(16) 



and expand the Bessel functions in equation (|T4[) in series 
in A. The resultant integrals can be evaluated analyti- 
cally knowing that 



da 



: K mWp 2 - Hpa) 

{p 2 - 2ipa) m / 2 



Then the implicit formula for A 



(17) 



3=0 



2pJ j\ K 2 {p) 



can be derived^. It can be rewritten as the explicit for- 
mula for A. The expression for the transverse permittiv- 
ity 

e T = l-A = l-fX 2( " +1) n%#7-^ (19) 



n=0 



2=0 



K 2 (p) {-2p)H\ 



can easily give the dependence K (CI). Taking n = one 
obtains 



2 K 2 , KijT- 1 ) 



(20) 



in a high-frequency limit O ^> 1, coincident with the 
expression in Ref. (10). As the frequency goes down, 
terms with higher n become important. The consecu- 
tive approximations are expected to work well for any 
temperature. The author believes this elegant derivation 
of the high-frequency approximation was not previously 
outlined. 



SL = 1 



2T(k- 1) 



t(k + i)r(« - §)*r 2 



r(«+l) (21a) 



0F(1- A)?- k A k - 1 T[k+1 



— kcsc(7tk) 



3 2 #i(l,™,2-s, A 
V^FA^CA-^^-W^+^j, 



et = 1 — 



r(«-i) 



2T(k+ |)r(*s- \)K 2 



2T(k + 1) (21b) 



20?A K (1- A)'- K r^/e+^ 



csc(7tk) 



+ ttk(2 - 3A - 2k) 2 Fi ( 1, --,2 — «, A 



2tt(Ak + k - l) 2 Fi ( 2, --,2 - k, A 



2V^A K (A- 1)^(71^+ i 



where u is determined by Landau rule (see equation ([8]) 
and discussion therein). The notation 2 F\ represents the 
regularized hypergeometric function. For the integer k 
the limit of k going to that singular integer value should 
be considered. The non-singular expressions were also 
derived, but are longer and are not provided here. The 
Mathcmatica 6 convention of branch cuts should be used 
to evaluate the values of roots and non-integer powers. 
It sets the branch cuts to be on arg(z) = n line in the 
complex z plane. 

The high-frequency limit for the transverse waves is 



tt 2 = K 2 



7rcsc(7rft)r(K) 



r(2 - k)t(k - § )r(« + i) 

Again the limit must be considered for integer k. 

2. Special cases 



(22) 



B. Power- law distribution 

The permittivity tensor for the power-law distribution 
(jTTJ) of particles can be calculated analytically. 



Much shorter formulas for permittivities can be derived 
for certain n. The shortest ones are for k in the observa- 
tionally motivated range k £ [2, 3]. Let me choose n = 2 
and k = 5/2 as the examples. 

The case n = 2 (equivalent to dNfa) ~ 7 -2 <i7 at high 



7 



7) leads to 

(4 - 10A)v / A^"T + 3ctttA 2 + 6A 2 arcsec(v / A) 



et = 1 • 



£L = 1- 



37rif 2 (A - l) 3 / 2 



(23a) 



+ 16A)VA - 1 + 12cr7rA 2 - 24Aarcsec(%/A) 
3vrif 2 (A - l) 3 /2 ' 

(23b) 

The case n = § (equivalent to dN(-f) ~ "f~ 3 dj at high 
7) leads to 

A(10-15A + 8A 3 / 2 (l-2a))-3 . 
^ + -- 16if 2 (A-l) 2 (24a) 



, E A(6 + 3A-8A 1 /2(i_ 2( j))- 1 
£i = 1 + 57r ^ 16/f 2 (A-l) 2 (24b) 



The longitudinal dispersion relation (|24bj) for k = 5/2 
remarkably gives the compact analytic form for fi(if). 
The undamped mode at low X exists in this case along 
with the single damped mode at higher K as 



(25a) 



for K = 0. 



/57T 



0= -y 48i4T 2 + 5tt - z 
for K = A^ ..I 



(16if 2 - 5tt) 



5tt 



(25b) 



4 ..-V5tt(9 + 4V5). 



The relation 17 (if) can be expanded near the point if 2 = 
V57T/4 to give the powers of (K x — if) 3 / 2 so that SIk 
stays real for K < K x and gains the imaginary part for 
K > K x . The second derivative of O(if) is discontinuous 
at K x . The explicit dispersion law O(if) ([25]) is plotted 
for k = 5/2 on Fig.QH 



C. Hybrid (3 distribution 

The analytic expressions for transverse and longitudi- 
nal permittivities need some terms to be abbreviated for 
brevity in the case of hybrid (3 distribution (fTB")) . They 
read 



e t 



1 



C{2AB + (3 4 ((B - 2C 2 ) (26a) 



ttA 3 BCK 2 ^ 

x arctan^l — i?arcsec/3)) — 2^4 3 A 2 arctanhC 

A 2 (3oi 
BCK 2 ' 




K 
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FIG. 12: (Color online) Longitudinal dispersion law for n = 
5/2: undamped mode (dot-dashed blue line), damped mode 
(solid red line), separated by Re(f2) = K (thin black line). 



£L = 1 



2/3 



B 2 C*/3 4 arcsec^-v4C (26b) 



nA 3 B 2 CK 2 
x (B(B + A A 2 ) + A/3 4 (B + 2C 2 )arcta,n A) 

3 A / o 1 r>2 o2i 



+ 2A 3 A( J B + C 2 /3 2 )arctanhC 

2f3Aai{B + C 2 /3 2 ) 
B 2 CK 2 



where 



A = ^J(3 2 - 1, B = /3 2 -A, C = \A~~A. (26c) 

The expressions for (3 = 1 coincide with those for the 
power-law distribution with k = 2. 



V. DISCUSSION & CONCLUSION 

The paper derives and reconsiders the broad range of 
linear wave effects in non-magnetized one species plasma. 
The various distributions are employed for comparison 
and to consider the realistic non-equilibrium plasmas. 
The thermal ([TO]) , power-law (|11[) . relativistic Lorentzian 
k (fT"2|) and hybrid (3 (fTB"]) distributions are employed with 
the possible inclusion of the simulated distributions^ in 
the future. 
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For longitudinal waves the maximum K was found to 
exist above that neither the undamped mode nor the 
damped mode survive. The mode completion effect 
was closed. The regions of the principal and the second 
damped modes overlap in the temperature-wavenumbcr 
plane for thermal plasma. The second damped mode al- 
ways exist. The lower possible K is zero for low tempera- 
tures for second and third damped modes, but is non-zero 
for higher temperatures. In contrast, for the power-law 
distribution (jlljl the second damped mode does not ex- 
ist for low k. When it exists at higher k, its lowest K is 
always zero. 

Some analytic results are derived that can accelerate 
the calculations and provide some insight. The full an- 
alytic result for longitudinal Q(K) was derived for the 
power-law distribution (fTTjl with n = 5/2. The transition 



at K x between the undamped mode and the principle 
damped mode was found to be smooth with discontinu- 
ous second derivative of Q(K) because of (K x — K) 3 / 2 
term. 
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